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摘要 : 为 了 探讨 我 国 不 同 地 区 菜 粉 蝶 种 群 之 间 的 遗传 分 化 情况 , 本 研究 运用 克隆 测序 法 ， 对 我 国 16 个 地 区 79 HK 
TE Pieris rapae rDNA 的 ITS-1 基因 序列 进行 测定 并 分 析 ; 同时 , DAZRZT STIR Pieris canidia HAME, 重建 了 它们 的 
系统 发 生 树 , 初步 探讨 了 它们 的 生物 地 理学 分 化 格局 。 序 列 分 析 结 果 显 示 : 所 测 菜 粉 蝶 的 ITS-1 序列 长 度 介 于 337 ~ 
458 bp 之 间 。 经 序列 比 对 后 , 在 347 个 位 点 中 共有 281 个 保守 位 点 , 63 个 变异 位 点 , 16 个 简约 信息 位 点 和 47 个 单 
突变 位 点 ; AMOVA 软件 分 析 其 核 苷 酸 多 态 性 指数 (nucleotide diversity) Pi 为 0.014, 每 位 点 Theta (per site) Eta 为 
0.041, 遗传 分 化 指数 (Fy) 为 0.351 (P «0.05) ; DnaSP 软件 共 检 测 出 48 个 单 倍 型 (haplotype) 序 列 , 单 倍 型 多 样 性 
( haplotype diversity) 的 频率 为 0.989。 系 统 发 生 分 析 结 果 表明 : 现存 的 菜 粉 蝶 各 地 理 种 群 与 其 地 理 分 布 没有 直接 的 对 
应 关系 根据 系统 发 生 树 结果 推测 , 该 蝶 种 在 我 国 的 最 早 建 群 可 能 发 生 于 辽宁 、 内 蒙古 及 北京 一 带 , 其 后 , 通过 幼虫 寄 
主 植物 的 水 陆运 输 以 及 成 虫 的 迁 飞 回 我 国 的 不 同 地 区 扩散 。 
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Abstract: In order to investigate the genetic variations of different geographical populations of Pieris rapae 
Linnaeus in China, the ribosomal DNA gene interval transcribed spacer 1 (ITS-1) sequences of 79 
indididuals from 16 areas were determined and analyzed. Upon these gene sequence data and using Pieris 
canidia as the outgroup, we reconstructed phylogenetic trees by means of neighbor-joining and minimum 
evolution analyses, aiming to clarify their biogeographical divergence pattern. The results of the sequence 
variation analyses showed that the sizes of the ITS-1 of Pieris rapae range from 337 to 458 nucleotides in 
length. Aligned with the ClustalW software, the total 347 sites, which contain 281 conserved, 63 variable 
and 16 parsimonious informative, were obtained. Meanwhile, for the whole population of this study, the 
nucleotide diversity Pi and the Theta ( per site) Eta were calculated to be O0. 014 and 0. 041, respectively. 
Moreover, total 48 haplotypes were identified and the haplotype diversity 0. 989 was obtained by the 
software DnaSP. Additionally, the fixation index ( Fr) 0.351 (P «0.05) was calculated by the analysis of 
molecular variance ( AMOVA ) software. The result of phylogenetic analysis indicated no pronounced 
correlation between their haplotypes and the geographical regions. Based on the phylogenetic trees, it is 


inferred that the Pieris rapae in China probably colonized first in the Liaoning, Inner Mongolia and Beijing, 
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and then spread to other areas by transportation of their host plants and migration of their adults. 


Key words: Pieris rapae; geographical population; genetic diversity; phylogeography; ITS-1 sequences 


3 4j HE Pieris rapae Linnaeus 隶属 于 昆虫 纲 
(Insecta ), 8$ WH H ( Lepidopetera ), 9j UE $} 
(Pieridae) , Ji E&Jg Pieris, 广泛 分 布 于 中 国 大 部 分 
地 区 。 其 幼虫 一 一 菜 青 虫 ,为 害 十 字 花 科 、 葡 科 、 旋 
AERE. EL ARRIR RRRA WES 9 科 35 种 植物 ( 李 
zw, 2006), 对 十 字 花 科 蔬 全 ,如 介 蓝 甘蓝 \ 油 荣 、 
花椰菜 的 危害 尤为 严重 。 目 前 该 种 蝶 类 的 研究 主要 
集中 于 形态 学 、 显 微 技 术 、 生 态 毒 理 \. 电 生理 研究 、 同 
工 酶 以 及 系统 发 生 等 方面 (Kim et al., 1988; 王 见 杨 
和 黄 可 威 , 2001; Cristian et al., 2002; 5H% E AIH 
哲 民 , 2002; ERE, 2005; Cerruti et al., 2007; 
许 丽 等 ,2007; 翟 兴 礼 和 杨 霞 ,2008; Jógar et al., 
2009) ,而 利用 分 子 遗 传 标 记 开 展 菜 粉 蝶 不 同 地 理 
种 群 遗 传 分 化 及 系统 地 理学 的 研究 鲜 有 报道 。 

利用 分 子 生物 学 技术 人 研究 不 同 地 区 采 粉 蝶 种 群 
的 遗传 学 关系 , 可 以 探索 各 地 菜 粉 蝶 种 群 在 数量 和 
空间 上 的 内 在 联系 ,从 而 为 制定 菜 粉 蝶 在 不 同 地 区 
的 区 域 性 控制 方法 及 沫 青虫 的 病虫害 防治 提供 科学 
依据 。ITS-1 序列 是 nrDNA 上 位 于 18S 5 5. 8S 
rDNA 之 间 的 序列 ,由 于 其 不 参与 核糖 体 的 形成 ， 
受 目 然 选 择 的 压力 较 小 , 进化 速率 较 快 ,因此 是 研 
究 种 间 和 种 内 分 子 进化 、 系 统 发 育 以 及 种 群 遗传 多 
样 性 分 析 的 最 常用 的 基因 之 一 (Schilthuizen et al., 
1995; ÆA, 2002; Annette and Victor, 2004) 。 

本 人 研究 尝试 运用 DNA 克隆 测序 法 对 我 国 16 个 
不 同 地 理 种 群 共 79 个 个 体 的 菜 粉 蝶 nrDNA ITS-1 
序列 进行 测定 和 比较 分 析 , 并 运用 系统 发 生 分 析 软 
件 重 建 其 系统 发 生 树 , 初步 探讨 该 种 蝶 类 在 我 国 的 
建 群 历史 及 其 随后 的 系统 地 理学 分 化 格局 ,以 期 为 
我 国 菜 粉 蝶 的 遗传 分 化 状况 和 今后 的 系统 地 理学 及 
相关 的 应 用 研究 提供 部 分 基础 资料 。 


1 材料 与 方法 


1.1 实验 材料 
材料 来 源 和 采集 时 间 见 表 1。 本 实验 所 用 标本 
经 鉴定 后 立即 放 和 人 9596 的 酒精 中 国定 , 市 回 实 验 室 
放 和 冰箱 中 - 207€ RF o 
1.2 基因 组 DNA 的 提取 、ITS-1 序列 的 扩 增 和 测定 
基因 组 DNA 采用 试剂 盒 DNeasy Blood & Tissue 
Kit (QIAGEN) 提取 ,所 用 引物 参考 于 等 (2003 ) 设 
计 的 1 对 引物 : ITS-1A: 5'-TACACACCGCCCGTCG 


CTACTA-3'; ITS-1B: 5'-ATGTGCGTICRAAATGTCGAT 
CTTCA-3'。 引 物 由 上 海 生 工 生物 工程 技术 公司 合成 。 

PCR 扩 增 体系 共 50 pL, 其 中 含 基 因 组 DNA 约 
0.5 ng, PCR 反应 程序 常规 PCR 产物 采用 
TRANSGEN 的 pEASY-T3 试剂 盒 进行 殉 隆 测序 , M 
序 工作 由 上 海 生 工 生物 工程 技术 公司 完成 。 
1.3 序列 处 理 和 系统 分 析 

用 BioEdit 7.0 ( Hall, 1999 ) 将 实验 克隆 后 测 出 
的 序列 进行 拼接 , 并 辅 以 手工 校正 。 在 GenBank 
( http://www. ncbi. nlm. nih. gov/) 上 确定 菜 粉 蝶 的 
ITS-1 序列 的 边界 及 长 度 。 用 Clustal X 1. 8 软件 
( Thompson et al., 1997) 进行 序列 比 对 。 核 酸 序列 
变异 及 单 倍 型 分 析 使 用 DnaSP 5. 0 (Librado and 
Rozas, 2009), LAZRZr SEX E Jg RE, Hj MEGA 4.0 
软件 (Tamura et al., 2007) 的 邻接 法 ( neighbor joining, 
NJ) 和 最 小 进化 法 (minimum evolution, ME) 重建 不 同 
单 倍 型 之 间 的 系统 发 生 树 。 采 用 Treeview 软件 
(Roderic and Edward, 1998) 显示 拓扑 结构 图 ， 利用 
Arlequin 3. 11 软件 (Schneider et al., 2007) 的 分 子 方差 
法 (AMOVA) 分 析 种 群 间 的 遗传 多 样 性 及 分 化 程度 。 


2 结果 与 分 析 


2.1 ITS-1 序列 组 成 及 其 遗传 多 样 性 参数 

本 人 研究 选择 齐齐哈尔 的 采 粉 蝶 和 西安 的 东方 菜 
粉 蝶 进 行 双 癌 克隆 测序 , 所 得 结果 用 BioEdit 7.0 进 
行 拼接 , 辅 以 手工 校正 。 并 在 GenBank 上 确定 其 
ITS-1 序列 的 边界 及 长 度 , 其 他 地 区 的 样品 与 之 对 
照 。 选 择 16 地 区 测序 较 好 的 79 条 ITS-1 序列 进行 
分 析 ( 表 2), 结果 表明 : ITS-1 序列 长 度 在 337 ~ 
458 bp 之 间 , 其 中 福建 三 明 ( 林 ) 菜 粉 蝶 ITS-1 序列 
最 短 , 辽宁 铁岭 (TL) 的 采 粉 蝶 序 列 最 长 ,而 外 群 辽 
TÀI (XC) 东方 菜 粉 蝶 仅 为 268 bp, 陕西 西安 
(XA) RIKERA 361 bp, 与 不 同 地 理 种 群 的 
羔 粉 蝶 有 了 明显 差异 , 其 中 北京 (BJ) 、 浙 江 椒 江 大 陈 
Ej ( DCD) 3 JH 9t BH ( 6Y ) , 181 JU f 2 2 PS ie PE 
( HEB) 安徽 芜湖 (WH) 和 山西 五 台山 (WTS) 6 个 
地 区 的 菜 粉 蝶 所 测 个 体 之 间 的 序列 长 度 完全 一 致 ， 
ITS-1 序列 呈现 较 高 的 稳定 性 ( 表 2), 而 辽宁 铁岭 
(了 TL) 的 菜 粉 蝶 序 列 间 长 度 差异 最 大 , 高 达 102 bp; 
除 此 之 外 ， 其 他 地 理 种 群 间 菜 粉 蝶 序 列 的 差异 很 
小 , 仅 为 1~4 bp。 序 列 经 Clustal X 1.8 多 重 比 对 
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R1 本 研究 使 用 的 莱 粉 蝶 地 理 种 群 样本 和 ITS-1 序列 信息 


Table1 Different Pieris rapae population samples and the GenBank accession numbers of their ITS-1 sequences in this study 


种 群 代码 采集 地 点 采集 时 间 "e "m 个 体 数 个 体 代码 CenBank 登录 号 
Population Sampling nn Number of Individual GenBank 
Collecting date Longitude Latitude 
code location individuals code accession no. 
WHI GU319925 
WH2 GU319931 
安徽 元 油 WH3 GU319924 
WH KE 2008.5 118. 42 31.33 6 
Wuhu, Anhui WH4 GU319931 
WHS GU319926 
WH6 GU319927 
BJl CU319889 
BJ2 CU319889 
北京 王家 园 BJ3 GU319889 
BJ REX 2009. 4 116.4 39. 88 6 
Wangjiayuan, Beijing BJ4 GU319889 
BJS GU319889 
BJ6 GU319888 
FJ1 GU319926 
FJ2 GU319900 
福建 三 日 FJ3 GU319901 
FJ me 5 2007.6 117.63 26.26 6 
Sanming, Fujian FJ4 CU319931 
FJ5 GU319931 
FJ6 GU319931 
CS1 GU319904 
TAE GS2 GU319935 
GS HRM 2007.6 106. 11 33.78 4 
Huixian , Gansu GS3 GU 3 19902 
GS4 GU319903 
GL1 GU319909 
GL2 GU319908 
西 
GL | WEM 2008. 8 110.29 25.25 5 GL3 GU319909 
Guilin, Guangxi 
GIA GU319907 
GLS GU319931 
GY1 GU319906 
GY2 GU319934 
贵州 贵阳 
GY 2008. 8 106. 66 26. 63 5 GY3 GU319905 
Guiyang, Guizhou 
GY4 GU319934 
GYS GU319934 
QQHEI GU319918 
黑龙 江 齐 齐 哈 尔 
QQHE n 2007.7 123. 82 47.35 3 QQHE2 GU319918 
Qigihaer, Heilongjiang 
QQHE3 GU319918 
HUB1 GU319913 
湖北 神农 架 ( 九 ; 
ARRE umh) HUB2 GU319934 
HUB Shennongjia ( Jiuchong) , 1998. 7 112.21 31.02 4 
Hubei HUB3 CU319913 
HUB4 CU319912 
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续 表 1 Table 1 continued 


种 群 代码 采集 由 ann ae uo MERO PERE 。 GenBank RE 





采集 时 间 经 度 纬度 
Population Sampling 集 7 Number of Individual GenBank 
Collecting date Longitude Latitude 
code location individuals code accession no. 
湖北 神农 架 ( 酒 过 坪 ) HEB1 CU319917 
HEB2 GU319915 
HEB Shennongjia ( Jiuhuping) , 1998. 7 112.21 31.02 4 
p HEB3 GU319916 
Hubei 
nd HEB4 GU319914 
DL1 GU319899 
辽宁 大 连 DL2 GU319934 
DL TAH 2008. 8 121.62 38.9 4 
Dalian, Liaoning DL3 GU319898 
DIA GU319897 
TL1 GU319922 
TL2 GU319922 
辽宁 铁岭 TL3 GU319921 
TL TW 2008. 7 123. 82 42.27 6 
Tieling, Liaoning TL A GU319920 
TL5 GU319923 
TL6 GU319919 
HHHTI1 GU319934 
HHHT2 GU319910 
蒙古 呼和浩特 HHHT3 GU319934 
HHHT HSCHTTRISHS 2008. 8 111.74 40. 83 6 
Huhehaote, Inner Mongolia HHHT4 GU319910 
HHHTS GU319911 
HHHT6 GU319910 
CF1 GU319891 
CF2 GU319890 
蒙古 赤 CF3 GU319892 
CF PIa Ri 2008. 8 118.9 42.24 6 
Chifeng, Inner Mongolia CF4 GU319893 
CF5 GU319894 
CF6 GU319910 
山西 五 台山 WTS1 GU319929 
WTS - 2008. 8 113. 32 38.72 3 WTS2 GU319930 
Wutaishan, Shanxi 
WTS3 GU319928 
YNI GU319931 
YN2 GU319932 
Z H YN3 GU319935 
YN 南 昆明 2009.4 102. 74 25.02 6 
Kunming, Yunnan YN4 GU319934 
YNS GU319933 
YN6 GU319934 
DCD1 GU319896 
ol DCD2 GU319934 
Wien 
DCD 2008. 7 121. 44 28.67 5 DCD3 GU319934 
Taizhou, Zhejiang 
DCD4 GU319895 
DCD5 GU319934 
XC1 GU319939 
洒 宁 兴 城 XC2 CU319937 
LL S 
XC 17 2008. 8 120. 73 40. 61 5 XC3 GU319938 
Xingcheng, Liaoning 
XC4 GU319940 
XCS GU319941 
陕西 西安 XAl GU319936 
XA ^ 2007.7 108.95 34.25 3 XA2 GU319936 





Xian, Shaanxi 


XA3 GU319936 
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后 ,发 现 部 分 地 理 种 群 内 和 种 群 间 序列 相同 ,其 中 
北京 (BJ) 浙江 大 陈 铝 (DCD) 38] JG T8 KA (HEB, 
HUB) 黑龙 江 齐 齐 哈 尔 (QQHE) 辽宁 铁岭 (TL)、 
广西 桂林 (GL) 的 种 群 内 部 分 序列 相同 ; 而 福建 三 
Bj CEJ) 贵州 贯 阳 (CGY) 内 花 古 呼和浩特 (HHHT) 、 
内 蒙古 赤峰 (CF) ZAE WH) HRE GS) 
等 地 理 种 群 间 部 分 个 体 序 列 却 存在 相同 性 ( 表 3) 。 
经 MEGA 4.0 软件 分 析 后 , 在 347 个 序列 比 对 位 点 
中 有 281 个 碱 基 序 列 为 保守 位 点 (conserved sites ) ， 
变异 位 点 (variable sites) 63 个 ,简约 信息 位 点 
( parsimonious informative sites) 16 个 和 47 个 单 突变 
位 点 (singletion sites) 。 经 DnaSP 5.0 软件 分 析 , 其 
核酸 多 态 性 指数 (nucleotide diversity) Pi 为 0. 014， 
每 位 点 Theta (per site) Eta 20. 041， 共 检测 出 48 
个 单 倍 型 ( haplotype ) 序列 ， 单 倍 型 多 样 性 
( haplotypes diversity) 频率 为 0. 989, 其 中 共享 单 倍 
型 共 3 种 : 其 中 DL2, YN6, YNA, DCD2, DCD3, 
DCD5, GY2, GY4, GY5, HHHTI , HHHT3 和 HUB2 
为 第 一 共享 单 倍 型 ( 记 作 Hap I); YNI, WE2, 
WH4, GL5, FJA, FJ5 和 FJ6 为 第 二 共享 单 倍 型 ( 记 作 
HaplI) ; QQHE1, QQHE2 和 QQHE3 为 第 三 共享 单 售 型 
( 记 作 Hap 焉 ) 。 而 独 享 单 倍 型 所 占 的 比率 以 及 分 布 频 
率 都 比较 高 , 各 种 群 独 至 单 倍 型 的 存在 及 分 布 表 明 这 
些 地 理 种 群 之 间 已 呈现 出 一 定 程度 的 遗传 分 化 。 

用 MEGA (version 4.0) 软 件 对 16 个 种 群 的 沫 
粉 蝶 的 碱 基 (T, C, A 和 C) 组 成 作 了 计算 , 发现 其 
碱 基 组 成 差异 不 大 , 79 个 个 体 的 T, C, A 和 G 的 平 
均 含 量 分 别 为 33.0% , 20.696 , 25.6% FI 20. 896 , 
G+C 的 含量 明显 低 于 A +T 的 含量 。 

表 2 莱 粉 蝶 与 东方 莱 粉 蝶 ITS-1 序列 长 度 
Table2 The nuclear ITS-1 sequence sizes of Pieris 
rapae Linnaeus and P. canidia Sparrman 


种 群 序列 长 度 (bp) 种 群 序列 长 度 (bp) 


Population Sequence length Population Sequence length 


BJ 340(6) HUB 339(3); 340(1) 

CF 339(1), 340(5) HHHT 339(2), 340(4) 

DCD 343(5) QQHE 2341(1),342(2) 

DL  382(2),342(2) TL — 356(3),357(1), 458(2) 
FJ  337(1),338(1),339(4) | WH 340(6) 

GL  342(2), 343(3) WTS  340(3) 


GS  338(1), 339(3) YN  337(1), 339(4) , 340(1) 
266 (1), 267 (1), 268 
(2), 272(1) 

XA — 356(1), 363(1) , 364(1) 

括号 内 数字 表示 序列 个 数 。Figures in brackets denote the number of 


sequences. 


GY  339(5) XC 





HEB  339(5) 
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RI 莱 粉 蝶 种 群 内 和 种 群 间 具有 相同 ITS-1 序列 的 个 体 列表 


Table3 The list of individuals with the same ITS-1 sequences 


within and between populations of Pieris rapae Linnaeus 


pe E 个 体 代码 
No. Number of Individual code 
individuals 

1 3 DCD2 , DCD3 , DCD5 

2 6 GY2,GY4,GY5, HHHTI , HHHT3 , HUB2 

3 2 QOHE2 , QOHE3 

4 3 YN1 ,WH2,WH4 

5 2 FJ4 , FJ6 

6 2 HUB1 , HUB3 

7 1 FJ1 , WHS 

8 2 GS2 , YN3 

9 2 TL1 , T12 

10 2 GL2 , GLA 

11 6 BJ1 , BJ2 , BJ3 , BJ4 , BJ5 , BJ6 

12 4 CF6, HHHT2 , HHHT4 , HHHT6 


2.2 种 群 间 ITS-1 序列 的 遗传 距离 

本 研究 采用 AMOVA 软件 对 菜 粉 蝶 的 单 倍 型 
( haplotype) 数据 进行 分 析 , 结果 显示 ,其 单 倍 型 多 
样 性 为 0 ~1.5 ZE, 平均 为 0.989, 其 中 浙江 大 陈 
(DCD) , 8 £& — B] (EJ) , x PM 5t EH ( GY) JER 
(BJ)、 内 蒙古 呼和浩特 (HHHT) 和 湖北 神农 架 九 冲 
(HUB ) 值 最 高 , 为 0.5 ~1.5; 而 黑龙 江 齐 齐 哈 尔 
(QQHE ) 最 低 ; 其 核 苷 酸 多 样 性 为 0 ~ 0. 96 之 间 ， 
平均 为 0.014, 其 中 湖北 神农 架 酒 索 坪 (HEB ) 值 最 
高 , 为 0.2 ~0.96 之 间 ，QQHE 最 低 ; 平均 核 甘酸 
差异 数 最 高 和 最 低 的 种 群 分 别 为 湖北 神农 架 酒 壶 坪 
(HEB) 和 黑龙 江 齐 齐 哈 尔 (QQHE) ( 表 4) 。 

结果 还 显示 : 菜 粉 屿 的 平均 遗传 分 化 距离 为 
0.351 (P «0.05) , 其 种 群 内 和 不 同 种 群 间 均 有 较 
高 的 遗传 变异 。 但 总 的 来 说 , 变异 主要 来 源 于 种 群 
内 ( 占 总 变异 的 64.93% ) ,其 种 群 间 和 种 群 内 的 遗 
传 分 化 指数 (Fsr) ILR 5. 
2.3 RANAR 

以 西安 (XA) 和 辽宁 兴 城 (XC) 的 东方 菜 粉 蝶 为 
外 群 , 分 别 以 邻接 法 (NJ 树 )、 最 小 进化 法 (ME ) 构 
建 了 16 个 地 理 种 群 共 48 个 不 同 单 倍 型 的 系统 进化 
树 , 重复 1 000 次 抽样 估算 系统 树 上 各 分 校 的 自 引 
导 值 。 结 果 显 示 两 种 树 的 拓扑 结构 基本 一 致 , 现 以 
NJ 树 为 例 说明 ( 图 1) : 各 地 理 种 群 均 具 有 和 较 高 的 单 
倍 型 多 样 性 , 有 旦 同一 种 群 的 不 同 单 倍 型 之 间 有 一 定 
HJ E, JÉ B AN fn] BJ M 2S DS AS. Horn ER 3 
(DCD) 桂林 (GL) , zx B8 (YN) 和 湖北 神农 架 九 冲 
( HUB ) 变异 最 为 明显 。 其 中 , 北京 (BJ) 呼和浩特 
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表 4 莱 粉 蝶 不 同 地 理 种 群 的 基因 单 倍 型 . 核 背 酸 多 样 性 .Tajima's D 中 性 检验 及 平均 核 音 酸 差 异 数 
Table4 The haplotype diversity, nucleotide diversity, Tajima's D neutrality test and mean number 
of pairwise nucleotide difference of Pieris rapae populations in this study 
FREE 单 倍 型 多 样 性 核 苷 酸 多 样 性 Tajima’s D 平均 核 苷 酸 差异 数 
Population Haplotype diversity ( Hd) Nucleotide diversity ( 7) Mean number of pairwise difference 
DL 1. 0000 x 0. 2722 0.015595 x 0. 012865 0. 00000 2.66667 
HEB 1. 0000 +0. 1768 0. 581952 x 0. 380973 1.43823 182. 33333 
WTS 1. 0000 x 0. 2722 0. 005900 x 0. 005562 0. 00000 2 
GS 1. 0000 x 0. 2722 0.011799 x 0. 010028 0. 00000 3. 33333 
YN 0. 7576 +0. 0927 0. 006818 x 0. 004510 — 0. 64365 1.65152 
DCD 1. 0000 x 0. 5000 0.011628 x 0.013000 0. 00000 4 
GL 1. 0000 x 0. 2722 0.019380 x 0. 015691 0. 00000 5. 33333 
GY 1. 0000 x 0. 5000 0. 005900 x 0. 007226 0. 00000 2 
FJ 1. 0000 x 0. 5000 0.014749 x 0. 016157 0. 00000 4 
TL 1. 0000 x 0. 1265 0. 152174 x 0. 093039 — 2.58268 7.2 
WH 1. 0000 +0. 1768 0. 010294 x 0. 007871 — 0. 06501 2.16667 
QOHE 0. 0000 x 0. 0000 0. 000000 x 0. 000000 0. 00000 0. 00000 
BJ 1. 0000 x 0. 5000 0.002941 x 0. 004159 0. 00000 1 
HHHT 1. 0000 x 0. 5000 0. 008824. x 0. 010189 0. 00000 3 
HUB 1. 0000 x 0. 5000 0. 023460 x 0. 024884 0. 00000 5 
CF 1. 0000 x 0. 1265 0.014118 x 0. 009687 —1.19955 4.4 
表 中 Hd 和 7 [J^ PE +R, The Hd and m values are mean + SD. 
3 5 不 同 地 理 种 群 间 莱 粉 蝶 的 遗传 分 化 指数 (Psr) 值 
Table 5 The F,r values of the 16 Pieris rapae populations in this study 
DL HEB WTS GS YN DCD GL GY TL WH  QQHE BJ HHHT HUB Fs 
DL 0. 505 
HEB 0. 257 —1.10 
WTS 0.450 0.254 0.528 
GS 0.391 0.251 0.182 0. 514 
YN 0.550 0.606 0.218 0.222 0. 520 
DCD 0.606 0.139 0.676 0.571 0.669 0. 521 
GL 0.486 0.254 0.435 0.385 0.522 0.345 0.496 
GY 0.352 0.116 0.143 0.032 20.105 0.625 0.327 0. 532 
FJ 0.410 0.119 20.357 0.250 0.434 0.571 0.384 0.222 0.516 
TL 0.356 0.391 0.339 0.337 0.614 0.230 0.291 0.231 0.245 — 0. 043 
WH 0.470 0.334 0.303 0.268 0.148 0.626 0.343 0.197 0.242 0.382 0.515 
QQHE 0.719 0.129 0.842 0.721 0.738 0.857 0.597 0.878 0.762 0.263 0.693 0. 542 
BJ 0.678 0.123 0.789 0.667 20.702 0.783 0.478 20.786 0.700 0.238 0.640 0.933 0. 537 
HHHT 0.619 0.123 0.687 0.587 20.663 0.696 0.398 0.667 20.600 0.230 0.567 0.800 0.333 0. 527 
HUB 0.389 0.112 0.333 0.250 0.336 0.455 0.103 0.167 0.235 0.197 0.042 0.500 0.182 0.000 0. 501 
CF 0.588 0.391 0.539 0.507 0.577 0.616 0.421 0.479 0.529 0.436 0.477 0.582 0.185 0.059 0.096 0.502 


最 后 一 列 为 种 群 内 Forto The Fs; values in the last column show those within each population. 


( HHHT) L TERIS (TL) 及 赤峰 (CF) 种 群 之 间 具 有 较 
近 的 亲缘 关系 , 与 其 实际 地 理 距 离 呈 现 出 一 定 的 相关 
性 , 但 是 这 些 种 群 内 的 部 分 序列 却 也 出 现 了 一 定 的 分 
化 , 如 CF5 5; YN1 和 WHI 具有 较 高 的 同 源 性 。 


3 讨论 


不 同 地 理 种 群 间 的 菜 粉 蝶 的 遗传 分 化 
群体 遗传 结构 是 指 遗 传 变异 在 物种 或 群体 中 的 


3.1 


一 种 非 随机 分 布 , 即 遗 传 变异 在 群体 内 、 群 体 间 的 分 
布 样式 以 及 在 时 间 上 的 变化 ( Hamrick and Loveless, 
1989) 。 这 种 非 随 机 分 布 是 由 不 同 过 程 共 同 作用 产生 
HJ, 包括 物种 长 期 的 进化 历史 、 分 布 区 改变 、 生 境 破 
碎 、 群 体 隔离 突变、 遗传 尝 变 繁育 系统 基因 流 与 选 
(Slatkin, 1987, 1995)。 因 此 , 要 阐明 生物 进化 
的 过 程式 样 和 机 理 , 必须 首先 探讨 生物 群体 遗传 变 
异 大 小 遗传 结构 及 其 变化 规律 以 及 影响 群体 遗传 结 


1150 昆虫 学 报 Acta Entomologica Sinica 53 3 


68 
46 


T] 


71 


72 
69 


80 


76 


74 


60 


99 98 
98 


98 


99 
45 


73 





71 GS4 


YN3, GS2 
GS1 

GY1 

DL4 
HEB4 
HEB2 
HAPI 


s4 [一 WIS3 


GS3 
DCD4 
GL2, GL4 
YNS 


67 DL3 


DL1 
GY3 
WTS2 
HEB3 


78 C WTSI 


HEB1 
WH6 


71 WHS, FJ1 


FJ2 

FJ3 

HUBI, HUB3 
WH3 

YN2 

HAPII 

WHI 

GL3 


78 [ HAPIII 


CF5 

DCDI 

TL5 

TL3 

TL1, TL2 

HHHT5 

BJ1, BJ2, BJ3, BJ4, BJ5 
HHHT2, HHHT4, HHH6, CF6 
GL5 

CFI 


84 


CF4 

85 CF3 
TL4 
HUB4 


BJ6 

TL6 

XA1, XA2, XA3 
XC3 

XCI 

XC4 


93 — XC5 


XC2 


图 1 3e 48 个 ITS-1 单 倍 型 的 NJ 系统 发 生 树 (以 东方 菜 粉 蝶 为 外 群 ) 
Fig. 1 Neighbor-joining tree of the 48 Pieris rapae ITS-1 haplotypes using P. canidia as the outgroup 


构 的 各 种 因素 (Hamrick，1983 ) , 

Wright (1978) 提出 遗传 分 化 指数 ( Fsr) 介 于 0 
~0.05 之 间 的 种 群 遗传 分 化 很 弱 , 介 于 0. 05 ~ 
0.15 之 间 的 种 群 遗传 分 化 中 等 , 介 于 0. 15 ~ 0. 25 
之 间 的 种 群 遗传 分 化 很 大 ,大 于 0. 25 表明 种 群 遗 
传 分 化 极 大 。 本 研究 表明 : 菜 粉 蝶 主 要 的 遗传 变异 
应 来 源 于 种 群 内 , 各 种 群 内 的 Psr 值 一 般 在 0.5 7c 
右 。 而 种 群 间 的 情况 比较 复杂 ， 其中, 北京 与 黑龙 
江 齐 齐 哈 尔 种 群 的 遗传 分 化 指数 (Fsr) 值 最 大 ,高 
150.933; 黑龙 江 的 齐齐哈尔 与 山西 五 台山 `\ 浙江 椒 


江 大 陈 咏 及 吐 州 贵阳 也 高 达 0.8 以 上 ; 而 湖北 神农 
架 酒 棍 坪 与 浙江 大 陈 岛 、 贵 州 喧 阳 、 福 建 三 明 、 黑 龙 
江 齐 齐 哈 尔 、 北 京 、 内 守 古 有 呼和浩特、 湖北 神农 架 九 
冲 ; 湖北 神农 架 九 冲 与 广西 桂林 、 内 聚 古 赤峰 和 呼 
和 浩特 ; 云南 昆明 与 贵州 贵阳 、 安 徽 范 湖 以 及 内 家 
十 亦 峰 与 呼 和 洛 特 的 sr 值 首 小 于 0.15; 另外 湖北 
神农 染 九 冲 与 内 聚 古 呼和浩特 的 Rsr 值 却 仅 为 0。 
除 此 之 外 , 辽宁 大 连 和 其 他 所 有 种 群 间 ; 湖北 神农 
架 酒 谈 坪 与 山西 五 台山 \ 甘肃 徽 县 云南 昆明 广西 
桂林 等 其 他 种 群 间 的 Fsr 却 大 于 0.25。 
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总 的 来 说 , 菜 粉 蝶 种 群 的 遗传 分 化 较为 明显 ， 
并 旦 主要 体现 于 种 群 内 。 其 各 地 理 种 群 间 遗传 差异 
形成 的 主要 原因 , 可 能 与 人 为 因素 导致 的 地 理 种 群 
间 基 因 流 的 有 无 和 大 小 有 关 。 同 时 , 本 人 研究 显示 ， 
海拔 高 度 是 限制 菜 粉 蝶 的 基因 交流 的 一 个 重要 因 
ZR. ABI AL TRA AER Ju Ih ANE SE PE Mb F HJ S fs [e] BJ 
HEC, Aper SHE BS SEE ERE T HH Sr B P 
作用 。 

物种 的 适应 能 力 、 生 存 能 力 和 进化 潜力 与 其 遗 
传 多 样 性 高 低 密切 相关 , 一 般 而 言 , 单 倍 型 间 的 核 
酸 多 态 性 指数 (nucleotide diversity) Pi 是 衡量 群体 
内 变异 程度 的 重要 指标 之 一 。 对 于 大 多 数 动物 来 
W, Pi 值 在 0.01 以 上 时 被 认为 变异 较 大 (Neigel et 
al., 1993) 。 本 实验 显示 ,， 荣 粉 蝶 单 倍 型 间 的 核酸 
多 态 性 指数 Pi H 0.014, 也 在 一 定 程度 上 说 明 其 各 
主要 种 群 间 的 遗传 多 样 性 较 高 。 

Pashley (1986) 通过 研究 表明 ,对 植 食性 昆虫 
来 说 , 异域 分 布 的 昆虫 产生 种 下 分 化 是 由 地 理 和 人 
为 的 双重 作用 决定 的 。 随 着 栽培 技术 的 发 展 和 交通 
的 便利 , 菜 粉 蝶 的 寄主 植物 之 间 的 交流 日 益 频 繁 ， 
而 菜 粉 蝶 也 随 着 寄主 植物 被 带 到 不 同 的 地 区 ,表现 
出 基因 型 相似 程度 升 高 。 同 时 人 研究 也 发 现 相同 地 域 
的 莱 粉 蝶 个 体 之 间 共 享 单 倍 型 出 现 的 机 率 相对 较 
高 ,而 跨 地 域 不 同 个 体 之 间 却 也 拥有 相同 的 基因 
型 , 单 售 型 的 出 现 没 有 一 定 的 规律 可 循 。 本 人 研究 的 
系统 树 拓 扑 结 构 显 示 , 不 同 种 群 间 的 系统 发 生 格 局 
与 它们 的 地 理 分 布 之 间 没 有 明显 的 关联 。 如 采 排 除 
人 为 因素 , 仪 徘 沫 粉 蝶 目 和 喘 的 迁移 很 难 形成 如 此 复 
杂 的 基因 流 。 因 而 ,人 为 因素 可 能 是 影响 且 粉 蝶 种 
群 遗传 结构 和 遗传 多 样 性 的 重要 因素 。 

SIS UA T 56 AERE Cruciferae) 为 主要 寄主 植 
9. KEMERIT Fr HH Brassica oleracea L. 2S5 
K, "HEI2ES UH HR DS mM EM TEREE, Z 
数学 者 认为 其 起 源 于 现存 的 一 种 一 年 生 杂 草 
Brassica oleracea var. sylvestris, 也 有 人 认为 其 起 源 
于 地 中 海地 区 4 种 野生 近 缘 种 的 复合 体 (王建 林 
等 , 2006) 。 甘 蓝 类 蔬 业 是 从 清朝 年 间 开 始 逐 步 传 
ARE, 在 我 国 的 栽培 历史 还 比较 短 ( 刘 瑞 和 王 超 ， 
2006) ， 而 荣 粉 蝶 的 另外 一 种 寄主 植物 一 一 大 日 沫 ， 
是 中 国土 生 土 长 的 一 种 蔬 某 。 有 研究 表明 : 美国 的 
KAIRIE A Gi TE 1620 年 以 后 随 十 字 花 科 植 物 由 欧 
洲 人 融入 的 ( 陈 学 新 , 1997) , JE F ESTE SER E tE, 
可 能 是 由 甘蓝 类 蔬菜 的 输入 而 进入 国内 的 。 根 据 图 
1 的 结果 推测 ， 沫 粉 蝶 在 中 国 的 建 群 事件 可 能 发 生 


TIR WREATH. BERE, 一 
部 分 回 北 蔓延 , 到 达 齐 齐 哈 尔 等 地 ; 一 部 分 向 西部 、 
西南 部 扩张 到 了 山西 甘肃 ` 贵 州 、 云南 等 地 ; 为 一 部 
分 沿海 延伸 到 沈阳 、 福 建 \ 浙 江 、 安 徽 广 西 等 地 。 
3.2 ITS-1 序列 在 蝶 类 系统 发 育 中 的 应 用 

KREN MEA RER, 其 病虫害 研究 在 
生物 学 , 尤其 在 分 子 系统 学 方面 的 研究 还 很 匮乏 。 
本 研究 结果 表明 , 不 同 地 理 种 群 间 的 菜 粉 蝶 rDNA 
的 ITS-1 序列 遗传 分 化 较 大 , 同时 存在 明显 的 种 群 
内 个 体 间 变 异 。 本 研究 表明 核 ITS 区 为 高 变 区 , 是 
探讨 种 内 种 群 间 遗传 分 化 和 系统 发 生 的 有 效 分 子 标 
记 , 同时 也 可 作为 相似 种 间 快 速 鉴定 的 分 子 标 记 
(本 研究 中 它 可 以 较 快 鉴定 形态 上 相似 、 较 易 混 清 
HYS) ER S RARR ) o 

本 人 研究 所 采用 的 试验 材料 虽然 涉及 的 地 理 范 围 
很 广 , 但 总 体 而 言 ,由 于 采样 方式 、 采 样 条 件 、 试 验 
设备 及 人 研究 时 间 等 诸多 因素 的 影响 ,可 能 会 对 茶 些 
试验 结果 产生 影响 。 为 了 更 全 面 苔 握 沫 粉 蝶 地 理 种 
群 空 间 分 布 规律 的 内 在 联系 , 有 必要 在 今后 的 研究 
中 进一步 增加 样品 的 数量 和 应 用 其 他 分 子 标记 ， 以 
期 进一步 佐证 本 人 研究 的 结果 ， JUS ER BU DX bs 
制 . 病 虫害 防治 及 分 子 系统 地 理学 研究 提供 更 为 全 
面 的 分 子 生物 学 数据 。 


致谢 ”感谢 齐齐哈尔 大 学 王 守 军 硕士 \ 西 北 师范 大 
学 交大 杰 教 授 等 在 采集 样本 过 程 中 提供 的 大 力 帮 助 
以 及 安徽 师范 大 学 项 贤 领 博士 在 数据 处 理 上 给 予 的 
悉心 指导 ; 中 国 科 学 院 动物 研究 所 武 春 生 研 究 员 帮 
助 鉴别 部 分 标本 ,并 在 成 文 过 程 中 给 予 宝贵 建议 。 
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